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The time signals relevant for nonresonant spectral hole burning, a pump-wait-probe method de- 
signed to investigate slow relaxation, are calculated. The step-response function following the 
application of a high amplitude ac field (pump) and an intermediate waiting period is shown to 
be the sum of the equilibrium integrated response and a modification due to the preparation via ac 
irradiation. Both components are calculated for a class of stochastic dipole reorientation models. 
The results indicate that the method can be used for a clearcut distinction of homogeneously and 
heterogeneously broadened susceptibilities as they occur in the slow primary relaxation of super- 
cooled liquids and other disordered materials. This is because only in the heterogeneous case is a 
frequency selective modification of the response possible. 

PACS Numbers: 64.70 Pf,05.40.+j,61.20.Lc 



Disordered materials such as glasses, spin-glasses, disordered crystals and proteins exhibit 
non-exponential relaxation behavior on macroscopic time scales pj. In the past decade sev- 
eral experimental techniques have been invented in order to investigate the detailed nature 
of the non-exponential primary response particularly of amorphous polymers and supercooled 
liquids Q, Q. These techniques allow to specifically select a (slow) sub-ensemble and after- 
wards monitor its relaxation. They have shown the existence of dynamic heterogeneities which, 
however, are not static but fluctuate in timeQ. Recently the technique of nonresonant spectral 
hole burning (NHB) has been developed and applied to study the primary relaxation of super- 
cooled liquids on the time scale of seconds |J. In the meantime NHB has also been used to 
investigate the slow relaxation of disordered crystals §, amorphous ion-conductors Q and spin 
glasses^. The interpretation of the results has mainly been guided by the fact that via the appli- 
cation of a large amplitude ac field of frequency O (10 _2 Hz- • • 10 2 Hz) the sample absorbs energy 
of an amount proportional to the imaginary part of the susceptibility evaluated at the pump 
frequency f2|J. It was assumed that a frequency selective modification of the spectrum should 
be feasible only if the response is given by a heterogeneous superposition of entities relaxing at 
different rates. However, the results of a recent numerical investigation of the application of NHB 
to a solvable glass- model [pX]] suggest that this view may have to be revised. Beyond phenomeno- 
logical interpretations a deeper unterstanding of a pump-wait-probe experiment like NHB can be 
expected from a thorough analytical analysis, which, however, has not been presented up to date. 

It is the purpose of this Letter to provide a sound response theory for NHB. The experimental 
protocol of NHB consists of the following procedure ||. First, a large amplitude ac pump field 
Ep sin (fir) is applied to a sample in thermal equilibrium for a time t p = n(27r/0),n = 1, 2, 3, • • •. 
After this pump period the system is left in zero field for a waiting time t w before the step 
response is monitored via application of a small amplitude dc field E$ at time t = 0. 

In contrast to ordinary nonlinear response theory, I calculate the linear response of a variable 
S(t) to the static field Eg but with a disturbed initial state of the system prior to the application 
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of Eg. This initial state is calculated in order 0(Ep). Correspondingly, the response reads as 
R*(Es, Ep,t,t w ,t p ) = R(Es,t) + AR(Es, Ep,t,t w ,t p ) where R(E$,t) denotes the equilibrium 
integrated response. The 'modification' AR(E$, Ep,t,t w ,t p ) originates from the deviations from 
thermal equilibrium at time t = and consists of two terms of order O(Ep) and O(Ep), respec- 
tively. A subtraction of the signals following the application of a positive and a negative dc field 
Es along with the addition of the signals obtained for a positive and a negative pump allow to 
extract the relevant terms in 0(Es),0(Ep). This corresponds to the phase cycle employed in 
experiments Q. The final result reads as: 

R*(t,t w ,t p ) = R(t) + AR(t,t w ,t p ) (1) 

with R(t) = ±(R(E s ,t) -R(-E s ,t)) = R(E s ,t) and 

Ai?(t,^,t p )=AA^rTr{s(0) e - j£ ('-^[-^ 1 ( J E 5 )]e- i£ ( T+ '^Ap 2 (t p )} (2) 

Here, C is the Liouvillian, C\{E$) the first order perturbation due to the dc field and M is a 
constant. Ap2(t p ) denotes the ©(i^p )-deviation from thermal equilibrium of the density matrix 
directly after the pump. This results from two terms, Ci(Ep) and £2(Ep) of respective order 
0{E P ) and O(Ep)^}. It is evident from Eq.@ that a non-vanishing AR(t,t w ,t p ) results only 
because the system has been driven out of equilibrium prior to the linear response experiment 
and Eq.(|l|) shows that the modified response is just the sum of the ordinary linear response and 
the modification. Therefore, in the NHB experiments the linear response starting from a modified 
inital state (in O(Ep)) is monitored as opposed to ordinary nonlinear (dielectric) experiments. 

In order to be able to further discuss the implications of Eqns.(|l],0), in the following I will 
consider explicit models of stochastic dynamics without inertial terms. This is a reasonable 
assumption for the slow relaxation processes of interest in the present context. The results 
obtained so far can be used directly for the corresponding Fokker-Planck (FP) or master equation 
(ME) obeyed by the conditional probability P(x, t), where x(i) denotes the stochastic process 



under consideration |12],||L3[|. For a FP equation, the nonlinear response theory can be formulated 
in general terms 0]. In this case the resulting expressions are particularly simple because higher 
order terms like £2(E) vanish. 

If a ME is considered instead the situation can be quite different. The procedure of calcu- 
lating the response as for a FP equation can be utilized only if the transition probabilities in 



the ME are chosen in such a way that a Kramers- Moyal expansion! 12] is possible. If, however, 
the 'jump length' in the transition probabilities is large, the question how to couple the tran- 
sition probabilities to the field arises. This can be understood most easily with the following 
argument. Denote the transition probability for a Xi — > x k transition without externally applied 
fields by Kk,i(t) = ^(xfc, a^; i), the time dependent probability distribution Pi(t) = P(xi,t) and 
the equilibrium probability by p eq = Pi(t — > oo) in a discrete notation. From ordinary statistical 
mechanics the change in the p eq oc e~@ e ' with (5 = \/{kpT) (kp is the Boltzmann constant) 
and €i the energy, is known. If the field E couples to a function M(x), these are just given by 
p e j q {E) = p^ q exp {(3M{xi)E). How to change the Kk,i(t) when an external field is applied is how- 
ever not easily determined in the general case. The only restriction is that detailed balance has 
to be obeyed and this will be used in all what follows. The most general form for the Ki^(E,t) 
therefore is 

Ki,k(E,t) = Ki,k(t)e? Ei5/i < k where * i)fc = aM( Xi ) - (1 - a)M(x k ) (3) 

with a denoting a real number. The perturbation series for the ME follows from Eq.(|2|) if one 
substitutes (—iC) by the master operator W and the perturbation {—iC n ) by a corresponding 
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y( n ) i n 0(E n ). The matrix elements of W and are given by: 

i 

(V (n) )i,fc = ^ [W)] B j«i,fc(t)*i,* " E «i,*(*)*i,fc| ( 4 ) 

The constant a is determined by calculating the long time limit of the pk(E,t) as these have 
to coincide with p, (E) in the presence of external fields. For example, for a ME that obeys a 



Kramers-Moyal expansion, one has a = 1/2[12]. 



To fix the notation, in the following I will consider the reorientations of rigid molecules with 
a permanent dipole moment \i. Thus S(t) in Eq.(||) is identified with /i(t), the stochastic process 
x(i) with the orientation w(t) = (<p(t),6(t),tjj(t)) in terms of Euler angles and M(wj) = Hi = 
//cos0j. For simplicity, correlations due to the dipolar coupling of distinct molecules and local 
field corrections are neglected throughout. A simple way to model molecular reorientations is to 
assume that they occur in a completely random way. This is usually termed as a 'random jump 
model' |fl5| . In this model the transition probabilites for a change in orientation from wi — ► Wk 
are independent of Wi, Wk and are given by Kk,% = T/N Vi, k = 1, • • ■ , N. (If reorientations on a 
sphere are considered, N = 8tt 2 .) The equilibrium probabilities are p^ q = 1/N and therefore the 
choice of a is not determined a priori. 

In the following, I will consider a slight modification of this random jump model, for which 
the corresponding ME can still easily be solved analytically. This model is defined by the choice 

« M (t) = Zp?T(t) with Z = J2p- q (5) 

i 

which depends only on the final orientation of the transition. From this fact it is evident that 
for this model a = 1. Additionally, a time dependent rate T(t) is considered because this will be 
used in the later discussion. A simple choice yielding a non-exponentially decaying response is 

r(t) = bt b ~ l /T b G [0, 1] (6) 

which reduces to a time independent rate r _1 for 6 = 1. Though the considered model is very 
simple minded, it still shows the characteristic features of any stochastic model. In models 
with more complex (time independent) K^i several rates A m will occur as the eigenvalues of the 
master operator W. The mean decay rate in such more sophisticated models can be defined by 
jfJ2 m ^m = jfY^i k K i,k- As this just equals r _1 for 6 = 1, cf. Eq.(||), the simple model on 
average also reveals the features of more realistic ones. Furthermore, if only two orientations are 
considered, N = 2, the time honoured double well potential (DWP) model with an asymmetry 
(ei — €2) is recovered. 

The response for the model defined by Eqns.(^J^) is now calculated according to Eqns.(|^j2|) 
using (||). The resulting normalized response is given by $(t) = e~* '* and 
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A$(t, t w ,t p ) = -±^-A(n, 6)(tVr)e-* ' T e- 1 ^ (7) 

with (A/i 2 ) = {jj?) — (/x) 2 , (fi n ) = X^PiV?- The magnitude of the modification at a given pump 
frequency, i.e. the 'excitation profile', is given by (t p = 2nir/Q) 

A(Sl,b) = (^} 2 (A^ 2 )(Vt)V'£ /t f tP dh ( 1 dt 2 sm(ttt 1 )t b f 1 sin (Ot 2 )4 _le ^ /T 

Jo Jo 

+\ (>>V> - (/^ 3 >) (6/r)e-*VW'Vsin(m 1 )¥r 1 e*V- ( 8 ) 
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Eq.(|8|) reveals one of the central results of the present calculations, namely the frequency selec- 
tivity of NHB. A(Q, b) vanishes for large and small burn frequencies 0. This directly implies that 
a non- vanishing modification of the response can only be achieved for f2 on the order of t _1 . The 
sign of A(Q, 6) is determined by the prefactor of the second term and therefore depends on the 
details of the assumptions for the p^ q . Note that for a vanishing (fi) = one has A(Q, b) = 0. 
Furthermore, for a DWP model, it is easily seen that ( / u) 2 (^ 2 ) = (/i)(/i 3 ) and the second terms of 
A(£l, b) vanishes. The same holds approximately, if only small asymmetries are assumed in the 
general case. Therefore, in the following calculations the second term of A(Q, b) will be neglected 
throughout. 

A transparent expression allowing a simple interpretation for the excitation profile is obtained 
for the special case 6=1: 

A(U, 1) ~ ^ /U ) 2 (A^) e "(Or)e"(2^r)(l - e~^ r ) (9) 

Here, e"(milr) = mf]r/(l + m 2 $7 2 r 2 ). Therefore, A(Q, 1) is determined by the absorbed energy 
and higher order terms, in accord with what is expected on physical grounds[||. 

In Fig.l the normalized excitation profile A(£l, b)/A max is plotted as a function of the pump 
frequency for 6 = 1 and 6 = 0.7 and one pump cycle. For 6 = 1, the excitation profile is 
considerably narrower than for 6 = 0.7 and also than the Debye susceptibility e"(f2T). This 
means that the selection is restricted to a narrow band of frequencies. Decreasing 6 gives rise to 
broader excitation profiles. 

Important features of the modification A<I>(i, t w , t p ) given in Eq.(^) are the following, (i) 
If A&(t, t w , t p ) < 0, the modified response <£* decays faster than the equilibrium response. In 
the approximation used in Eq.(|9|) this always holds. The fact that A<$>(t,t w ,t p ) ^ only for £2 
on the order of r _1 has already been discussed above, (ii) AQ(t, t w , t p ) is only non-zero in a 
finite intervall of time t, determined by 6. The regime of A<i?(i, t w , t p ) ^ becomes larger with 
decreasing 6 (i.e. the 'spectral holes' become broader). The maximum value of the modification 
is found at the time t max = (l/r) 1 / 6 . For a homogeneous response (a fixed value of r) a variation 
of the pump frequency will only alter the overall amplitude of the modification. Therefore, it 
is not possible to perform a frequency selective modification in this case, (iii) The modification 
dereases as a function of the waiting time t w . This re-equilibration proceeds with the same 
relaxation time r. The reason for this behavior becomes evident from the form of the modified 
orientational distribution after the pump, p]f od (t p ) ~ (fik ~ {p))p £ k \ which in turn relaxes with 
the dipole relaxation time r. For more sophisticated models, other relaxation times A" 1 may 
occur here. The important finding is that there is no extra time scale for re-equilibration. It 
should be mentioned that similar results are found for a variety of stochastic models of dipole- 



reorientation[16] 



In order to further clarify the implications of the calculations, in the following I will consider 
specific examples. Often a stretched exponential function is used to parametrize the equilibrium 
response, <3?(t) = e~ <yt l TK ^ iK with (3k smaller than unity, see e.g. (p. In a heterogeneous scenario, 
this response is viewed as originating from a distribution of t§, i.e. = / dr 5 (r)e"*/ T with 
some distribution function g(r), whereas $(i) is assumed to be intrinsically non-exponential in a 
homogeneous scenario. It is instructive to compare directly the results of calculations for the two 
extreme scenarios for a realistically broadened susceptibility. Such a comparison of A$(t) is shown 
in Fig. 2 for (a) a heterogeneous and (b) a homogeneous scenario with <I>(t) = exp (— (t/ls) ). 
This means, I used 6=1 and a corresponding distribution of r in the calculations for case (a) 
and 6 = 0.7, T~ b = Is in (b). In the heterogeneous case it is evident that the time at which the 
maximum modification shows up strongly depends on the burn frequency f2. This demonstrates 
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the frequency selectivity of NHB. Such a dependence is missing completely in the homogeneous 
case. In the heterogeneous case the maximum modification shifts towards longer times with 
decreasing J7. This is because those relaxation processes with a r yielding the maximum A(Q, 1) 
(cf. Eq.(||) and Fig.l) contribute most to A <]>(£). The functional dependence of A$> max = 
maxitw — 0>ip) upon f2 depends sensitively on the choice of the relaxation time distribution 
(e.g. on the value of (3k)- Generally, there will be a rather strong dependence for large 0, which 
crosses over to a Sl-independence of A$ max when reaches the smallest relaxation time of the 
distribution. Exchange processes, which may be responsible for fluctuations of relaxation rates, 
may partially suppress the O-dependence . 

Also included Fig. 2(a) are experimental data on propylene carbonate^. It is remarkable that 
all main features of the data can be described by the simple heterogeneous model. The position of 
the modifications and the relative amplitudes are in quantitative agreement with the data. Also 
the width and the asymmetry are described with high accuracy. A similar quantitative agreement 
is obtained for the t^-dependence of experimental data on supercooled liquids [jO]]. 

So far, I have considered stochastic models for dipole reorientations. All these models show 
qualitatively the same features. In particular it is found that the re-equilibration during the 
waiting time t w takes place on the time scale of the intrinsic relaxation time (r in the above 
examples). In an experiment on a relaxor ferroelectric a re-equilibration time much longer than 
has been found 0. In order to investigate the theoretical conditions for such an effect I 
consider a simple 'nonlinear Debye model'. A physical realization of this scenario may be a 



system where the dynamics is dominated by domain wall depinning[17|. The polarization is 
assumed to relax with a correlation time To = T^e^ 11 where U a is the activation (pinning) energy 
and 1 / Too an attempt frequency. Application of an electric field will change the activation energy 
roughly by an amount U a (t) = U a (0) — cP sp E(t), where c is a constant and P sp , the spontaneous 
polarization^], which in the simplest case is proportional the the applied field, P sp , oc E(t). The 
relaxation of the deviation from thermal equilibrium of U a (t) is assumed to obey a linear law: 
dtSU a (t) +'ySU a (t) = cdt [E(t) 2 \ (meaningful for 7 < 1/tq), where the rate 7 may depend on the 
magnitude of deviation from equilibrium. This expression together with the equation for P(t), 
d t P(t) + (l/r)P(t) = x(^/r)E(t), can be solved in 0(E S ) and 0(E%) for the NHB situation. 
Here, % denotes the dielectric susceptibility and r = Toe SUa( - t \ The calculation leads to results 
similar to those obtained above with &(t) = e~ l ^ T ° and A$>(t,t w ,t p ) oc e~" ,tw . The important 
point is that in this model I/7, and not tq is the relevant time scale for re-equilibration. 

In conclusion, I have shown that NHB can be understood as a linear response experiment 
starting from a nonlinearly perturbed initial state. Generally, it is always possible to separate 
the linear response from the effects of the pump process. For systems with stochastic dynamics 
the method is clearly able to discriminate between homogenously and heterogeneously broadened 
dielectric spectra. A frequency selective modification of the response is possible only in the 
heterogeneous case. The re-equilibration during the waiting time is determined by the relaxation 
of the modified orientational distribution created during the pump period. The hole recovery 
does not show the appearence of a second time scale. A longer time scale for re-equilibration, if 
it is observed, has to be attributed to intrinsic non-equilibrium effects. 

I am grateful to R. Bohmer, H. Sillescu, O. Kircher, G. Hinze and R. Schilling for very fruitful 
discussions and comments. 
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Figure captions 

Fig.l : Normalized 'excitation profile' A n (tt,b) = A(£l,b)/A max versus £l/£l max for b = 
l,tt max T = 0.736 (full line) and b = 0.7,0 

max T — 0.185 (dashed line). Also shown is 
the imaginary part of the susceptibility, e"(fir), where r is the relaxation time, (dotted 
line) for comparison. The lines for A n (£l,0.7) and e"(Qr) have been shifted by 0.5 and 1.0 
units, respectively. 



6 



Fig. 2 : 10 3 A3>(i) versus rescaled time t/t max for (a) a heterogeneous and (b) a homogeneous sce- 
nario for various burn frequencies $1 and j3 2 Ep / (Afj?) = 0.05. In both cases the equilibrium 
response decays as $(t) = e _ ^ ls ^°' 7 . The time t max = 1.0s in case (a) and t max = 0.095s in 
case (b). The used frequencies are: heterogeneous scenario: ftr = 5.0 (1), 1.0 (2), 0.2 (3) 
0.1 (4); homogeneous scenario: $7r = 5.0 (1), 1.0 (2), 0.1 (3), 0.05 (4); Experimental data 
in (a) are adapted from ref.|§. Here, Or=1.02 (2) and Or=0.203 (3). 
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